arXiv:1501.05648vl [astro-ph.HE] 22 Jan 2015 


PROCEEDINGS 

0F SCIENCE 



Gamma-ray lines from SN2014J 


Thomas Siegert*and Roland Diehl 

Max Planck Institutfur extraterrestrische Physik, D-85748 Garching, Germany 

E-mail: tsiegert@mpe.mpg.de, rod@mpe.mpg.de 

On 21 January 2014, SN2014J was discovered in M82 and found to be the closest type la super¬ 
nova (SN la) in the last four decades. INTEGRAL observed SN2014J from the end of January 
until late June for a total exposure time of about 7 Ms. SNe la light curves are understood to be 
powered by the radioactive decay of iron peak elements of which 56 Ni is dominantly synthesized 
during the thermonuclear disruption of a CO white dwarf (WD). The measurement of y-ray lines 
from the decay chain 56 Ni—» 56 Co—> 56 Fe provides unique information about the explosion in su¬ 
pernovae. Canonical models assume 56 Ni buried deeply in the supernova cloud, absorbing most 
of the early y-rays, and only the consecutive decay of 56 Co should become directly observable 
through the overlaying material several weeks after the explosion when the supernova envelope 
dilutes as it expands. 

Surprisingly, with the spectrometer on INTEGRAL, SPI, we detected 56 Ni y-ray lines at 158 and 
812 keV at early times with flux levels corresponding to roughly 10% of the total expected amount 
of 56 Ni, and at relatively small velocities. This implies some mechanism to create a major amout 
of 56 Ni at the outskirts, and at the same time to break the spherical symmetry of the supernova. 
One plausible explanation would be a belt accreted from a He companion star, exploding, and 
triggering the explosion of the white dwarf. 

The full set of observations of SN2014J show 56 Co y-ray lines at 847 and 1238 keV, and we 
determine for the first time a SN la y-ray light curve. The irregular appearance of these y-ray 
lines allows deeper insights about the explosion morphology from its temporal evolution and pro¬ 
vides additional evidence for an asymmetric explosion, from our high-resolution spectroscopy 
and comparisons with recent models. 
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1. Introduction 

Supernovae of type la are thought to be standard candles for the cosmic distance ladder at dis¬ 
tances of more than sUOO Mpc. This is based on the uniformity of most measured light curves in 
the optical bands [1], It is not clear what the progenitor system of a SN la is, except that a CO white 
dwarf and its thermonuclear disruption must be part of it. There is more than one possible branch 
for a WD to enter this final fate. One possibility is a giant star companion delivering the additional 
fuel for the WD to ignite in its center when reaching its Chandrasekhar limiting mass [2, 3]. It can 
also undergo a triggered explosion by an accreted He layer which is detonating on the surface [4], 
These scenarios are called single degenerate (SD) and have been the consensus for several years. 
Recently, additional observations and increasing model sophistication revealed more diversity and 
offered plausible alternatives. A broader range of binary systems and more methods of igniting a 
white dwarf, independent of its mass, manage to destabilize the degenerate matter, e.g. by accretion 
flow instabilities [5] or by merging two white dwarves (double degenerate, DD) [4, 6, 7, 8]. The 
common product of all the different models is the radioactive isotope 56 Ni which is produced in 
nuclear fusion reactions at high densities, once the WD is ignited [9]. This explosive nuclear burn¬ 
ing competes with the expansion of the star, resulting in some outer parts of the WD being burnt 
towards intermediate mass elements instead of to iron group nuclei, or even left unburnt as carbon 
and oxygen [10]. Depending on the models and by indirect estimations [11], the typical amout of 
56 Ni produced is expected to be between 0.1 and 0.9 M 0 embedded in about 0.5 to 0.9 M 0 of other 
material [10, 12]. 

Up to several weeks, the SN cloud should be opaque even to penetrating y-rays from the 56 Ni 
decay with a characeteristic life-time of 8.8 days [13]. During this time, the outer gas absorbs the 
radioactive energy input and re-radiates it at lower energies from UV through IR, with such absorp¬ 
tion becoming less efficient as the SN expands. Neither the explosion dynamics nor the evolution 
towards explosions from WD properties and from their interactions with their companion stars can 
easily be assessed by observations in these wavelengths [14]. Only pre-explosion data brought 
some insights towards the nature of the binary companion star [3, 15, 16, 17, 18], and from its 
interactions with the SN [19, 20, 21], 

Deeper insights can be obtained by direct measurements of y-rays from the 56 Ni decay chain 
56 Ni—» 56 Co—)- 56 Fe, leaking out of the expanding SN after several weeks. According to most re¬ 
cent models, the consecutive decay from 56 Co to 56 Fe with a characteristic time scale of 111 days 
should show a maximum in y-ray emission brightness at days 60-100 after the explosion, declining 
thereafter due to the radioactive decay of 56 Co [22, 23, 24, 25, 2, 10, 26, 27, 4, 28, 29]. The bright¬ 
ness of 56 Co decay y-rays during the time after the explosion is determined by the initial amount 
of available 56 Ni and how it is distributed within the expanding SN. This is the phase in which the 
rise and fall of the y-ray line intensity provides unique information on the type of the explosion and 
the structure of the SN and thus gives hints on the progenitor system. 

SN2014J is a unique opportunity to investigate these aspects using the spectrometer SPI [30] on 
INTEGRAL [31] because it is the closest SN la in the last 40 years [32]. It was accidentally dis¬ 
covered by an undergraduate class lead by Steve Fossey on January 22, 2014 [33], and recognised 
as a type la explosion from early spectra [32]. It occured in the nearby starburst galaxy M82 at a 
distance of «3.3 Mpc [34], and is seen only through major interstellar material in the host galaxy 
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along the line-of-sight, causing rather large reddening [21]. This naturally involves difficulties in 
measuring the proprietary SN brightness and low-energy spectra at high precision. 

INTEGRAL started to observe SN2014J on January 31 [35], and pointed to it until June 26 ac¬ 
cumulating a total exposure of «7 Ms. This corresponds to observations of SN2014J during days 
16.3 to 164.0 after the inferred explosion date of January 14, UT 14.75 ± 0.21 [36]. 

Using our new analysis method for determination of the dominating instrumental background [37, 
38], for the first time we find SN la y-ray lines at 158 and 812 keV originating from the decay of 
56 Ni. At only 20 days after the explosion, the detection itself and the spectral shapes are unexpected 
and uncover asymmetries in the explosion. INTEGRAL’S follow-up observations have been used 
to provide a SN la 56 Co y-ray light curve. This, at first glance, conforms to most models. But with 
SPI’s unique spectral resolution (~2 keV), the observed 56 Co y-ray line shapes provide evidence 
for an asymmetry in SN2014J and also supports a double detonation scenario. 

2. A new analysis method 

Generally, in our spectroscopic analysis we fit scaling factors of proposed sky model inten¬ 
sity distributions plus models of instrumental background (BG) components to the set of spectra 
accumulated during the observations. This resembles linear regression of specific, given time se¬ 
quences, incorporating the instrumental (imaging) response function. Hence, data dy of energy bin 
k are modeled as a linear combination of the Nj sky model components M\j, to which the instrument 
response matrix R^ is applied, and the Nb BG components B /k : 

Ni Ni+N b 

dk = 'E R ik'L0i M U+ L Mjk (2-1) 

j i=l i=Ni+l 

i.e. the comparison is performed in data space with respect to time 1 , which consists of the counts 
per energy bin k measured in each of SPIN detectors j for each single telescope pointing as part of 
the complete observation. 

For several years, the common approach was to treat the count rate in each energy bin separately, 
thus not further using spectral response information of the instrument. The BG was determined 
empirically by using other instrumental rates as a kind of radiation monitor for the BG stimulating 
events. Since these rates do not exactly resemble the BG count rate over a longer period of time, 
additional scaling parameters had to be introduced. Moreover, the BG shows strong variations with 
energy, and for this reason these methods need to be calibrated for each energy bin separately. 

Our new approach for determining the BG count rate is based on detailed spectral fits in the energy 
region of interest. We thus investigate the nuclear reaction physics happening inside the satellite 
by studying the raw spectra for each detector and each orbit that is part of the analysed data set. A 
three-day integration is used to avoid statistical limitations in the BG count prediction per pointing 
and helps to identify even minor spectral features. We confirmed independently that spectral BG 
features do generally not change over such periods. We determine spectral shape paramters of each 
instrumental feature (see figure 1) for all detectors by fitting the spectra using a Metropolis-Hastings 

’The term "time" is somewhat misleading since one and the same pointing can be observed at different times. 
However, for simplicity, "time" is used in order to describe any means of "a particular pointing at time t". 
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Figure 1: Raw orbit-integrated count spectrum of the SPI camera. Shown are the photon counts in the energy 
band 780 - 920 keV per half keV energy bin in black and superimposed a spectral fit in red, identifying 28 
different components: A y-ray continuum, described by a power-law, and 27 instrumental y-ray lines (20 
already mentioned in e.g. [39], indicated with arrows). 

algorithm in a Monte-Carlo Markov-Chain and fix these spectral templates during one orbit. The 
spectral response function used for fitting the data is represented by a symmetric Gaussian func¬ 
tion, convolved with an exponential tail function towards lower energies which incorporates the 
degradation due to cosmic bombardment in each detector. It has four fit parameters per detector 
and orbit in each y-ray line [40]. Since each detector shows a slightly different degradation at a 
particular point in time, and in addition reacts differently to a certain BG source inside the satellite 
(being more or less illuminated than another detector), a unique and independent time-sequence 
of BG counts can be predicted. Long-term investigations showed that the relative intensity of a 
particular spectral BG feature (e.g. a y-ray line originating from the anti-coincidence shield) in one 
detector with respect to the mean intensity of the same feature in all detectors (detector ratios) is 
constant in time, and thus maximally independent of possible sky signals. This spectral-shape and 
detector-ratio constancy is imprinted onto the short-term variations while the mean amplitude (over 
all detectors) is still allowed to vary. 

Together with plausible candidate celestial source models, in the case of SN2014J one single point 
source, the scaling parameters 0, in equation (2.1) for both, BG and sky, are then determined by a 
maximum likehood fit simultaneously, setting up the BG and celestial spectra. 


3 56pqj y. ra y measurements 

Figure 2 shows the resulting spectra of a point source located at the position of SN2014J (at 
galactic coordinates (l/b) = (141,427°/40.558°)) at days 16.6 to 19.2 after the explosion. The 
two major lines at 158 and 812 keV can well be described as Gaussians on top of a constant 
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Figure 2: Spectra (black crosses) from the position of SN2014J between days 16.6 and 19.2 around 158 
(left) and 812 keV (right), showing the two main 56 Ni decay y-ray lines. The solid black lines are Gaussians 
as fitted to the data points. The scaled raw spectrum is superimposed in gray, together with the peak positions 
of the fitted Gaussians compared to the peak positions of the nearest instrumental lines. 

offset. Within the uncertainties, the two lines are identical. Also the imaging analysis locates the 
signal at the position of SN2014J, within spatial-resolution uncertainties, as shown in figure 3. 
Integrating the Gaussian line profiles yields 
the y-ray line fluxes to be (1.10 ± 0.42) • 

10~ 4 ph cm~ 2 s _1 for the 158 keV line and 
(1.90 ±0.66) • 10~ 4 ph cm -2 s” 1 for the 812 keV 
line. The fitted centroids for the two lines are 
159.43 ±0.43 keV and 811.84 ±0.42 keV, re¬ 
spectively, consistent with the laboratory ener¬ 
gies of 158.38 keV and 811.85 keV, thus corre¬ 
sponding to a bulk velocity below »2000 km s _1 
(2a conf.). The measured intrinsic celes¬ 
tial line widthssuggest no significant Doppler- 
broadening. Although broader components 
could underlie the main signals, the velocity 
spreads are below 2000 km s' 1 (2a conf.). The 
measured flux values can be converted to a 56 Ni mass seen to decay at ~ 18 days past explosion 
that corresponds to (0.06 ± 0.03) M 0 of 56 Ni at the time of the explosion. 

These values are rather small for typical scenarios in which 56 Ni is buried deeply in the SN cloud 
where the all 56 Ni y-rays should have been absorbed at this time. Our measurement implies that a 
substantial fraction of 56 Ni must have been produced close to the surface of the WD, at a depth not 
exceeding a few g cm~ 2 in column density. Such an amount of 56 Ni freely exposed, together with 
the constraints on kinematics, is unexpected in most explosion models. 

Observations in other wavelengths also suggest SN2014J to be an unusual explosion, since the ex¬ 
pected t 2 -behaviour [41] in the rise of the optical lightcurve is not strictly represented by the data: 
The light curve hours after the explosion is rather well described by a broken powerlaw, showing a 
steeper rise and a possible shoulder which might also imply 56 Ni near the surface [21, 36]. 
Classical SD scenarios have been excluded from upper limits of X-ray emission [42]. Sub- 
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Figure 3: Location of the 56 Ni line signal on 
the sky. Within instrumental uncertainties (la 
dashed, 2a solid circle), is the y-ray emission 
from the 158 and 812 keV line mapped onto the 
position of SN2014J (cross). 
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Figure 4: Sketch of an ejecta configuration compatible with our observations. Helium accreted in a belt 
before the explosion produces the 56 Ni belt at the surface of the ejecta. The 7 -rays can escape from the belt 
material, while the 56 Ni in the core (black) is still buried at high optical depths. The observer is believed to 
view the SN nearly face-on. 

Chandrasekhar models with a He donating companion star or DD merger seem to work better 
in the case of SN2014J and already explain the 56 Ni located near the surface [43]. Most classical 
double-detonation models predict high velocity 7 -ray lines to be observed if a shell of 56 Ni was sur¬ 
rounding the SN ejecta [44, 45, 7], Moreover, this configuration should have an imprint on optical 
observables [46] which are not seen in SN2014J [42]. A plausible, though somewhat speculative 
configuration that explains the 7 -ray data, is an equatorial He belt rather than a shell accreted from 
a He donor. If this situation would be observed essentially pole-on, our kinematic constraints can 
be met. Analogous to classical novae, where this idea has already been discussed [47, 48], unstable 
mass transfer on the Kelvin-Helmholtz time-scale from the He companion to the WD is possible 
if the He donor is more massive than the WD and fills its Roche lobe. In order to conserve a belt 
and not to spread the accreted He over the WD before the He flash explosion, either the accretion 
rate must be of the order 10 5 M 0 yr 1 [49] and/or the WD must be rapidly rotating so that the 
accretion is faster than it loses angular momentum. This high accretion rate can be achieved if 
the material is mostly He and not H. This would be the standard scenario for a supersoft X-ray 
source and is excluded for SN2014J [50, 42]. When this He belt becomes dense enough, explosive 
He burning may be ignited, and the configuration as shown in figure 4 can be established. Re¬ 
markably, this configuration is consistent with both, optical and 7 -ray observation data. Moreover, 
our radiation transfer simulations [37] do not produce easily distinguishable features in UVOIR 
wavelengths from the 56 Ni-belt compared to standard models which makes SN2014J appear quite 
normal, at any observations angle. 

4. 56 Co 7 -ray measurements 

Figure 5 shows the resulting spectra of a point source located at the position of SN2014J 
between days 66.3 and 99.1. Two 56 Co 7 -ray lines are clearly detected, and the flux ratio de- 
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Figure 5: Spectra (black crosses) from the position of SN2014J at the expected 7 -ray line intensity max¬ 
imum between days «60 and «100 around 847 (left) and 1238 keV (right), showing the two main 56 Co 
decay 7 -ray lines. The solid black lines represent Gaussians as fitted to the dervied data points. Using 
the information from this broad binning, the lines appear broadened, equivalent to a velocity spread be¬ 
tween 3000 and 8000 km s -1 , and blue-shifted about 3000 km s . The measured line intensities are 
(3.65 ±1.21) • 10 ~ 4 ph cm ” 2 s~ ! (847 keV) and (2.27 ±0.69) • 10 ~ 4 ph cm “ 2 s ” 1 (1238 keV), respectively, 
which corresponds to an initially synthesized 56 Ni mass of (0.50±0.12) Mq, assuming the SN is transparent. 

rived from the line intensities of (3.65 ± 1.21) ■ 10 4 ph cm ~ 2 s ” 1 from the 847 keV line and 
(2.27 ±0.69) • 10 ~ 4 ph cnT 2 s ” 1 from the 1238 keV line is 0.62 ±0.28 (lab value: 0.68). The mass 
of 56 Ni derived from these flux values is (0.50 ±0.12) M 0 , assuming the SN is fully transparent. 
This mass estimation, the constraints on the Doppler-broadening of about 40 to 50 keV (FWHM), 
corresponding to velocities of (5500 ±2500) km s' 1 , and the constraints on the position of the 
lines to be blue-shifted about (3000 ± 1500) km s -1 are in good agreement with the majority of SN 
la model calculations. 

However, our newly developed analysis method allows us to investigate the spectral shapes of the 
56 Co lines in even more detail. Figure 6 shows the same spectrum shown in figure 5 (left panel) 
at different energy binnings. Apparently, a single broad Gaussian does not capture the line shape 
properties any more, suggesting SN2014J to be an asymmetric explosion event. 

In general, the 56 Co 7 -rays should emerge gradually from the inner part of the SN as its envelope 
gets more and more transparent over time while expanding [51]. Thus, the spectral shapes at early 
times with respect to late times may be different, e.g. moving from a general blue-shift (near side) 
towards the laboratory energy. If ejecta would include major non-sphericities, velocity signatures 
could be more complex. Therefore, we separated our observations for both lines in four different 
epochs as shown in figure 7: The first epoch, ranging from days 16.3 to 41.3 after the explosion, 
includes the optical light curve’s maximum in which most 7 -rays from the inner core should be 
absorbed. Epoch (2), between days 41.3 and 66.3 is the phase in which the 7 -ray emission should 
gradually rise and end up at the maximum in the third epoch from days 66.3 to 99.1. After an 
observation gap between days 100 and 134, the fourth epoch (days 134.8 to 164.0) should capture 
the 7 -ray emission in a rather transparent SN. The epochs are chosen to be sensitive to different SN 
la model scenarios, to help distinguish among them. Epoch (4) is almost independent of the early 
evolution and traces the global dynamics of the 56 Co produced in the explosion. It should follow 
radioactive decay, in particular after the 7 -ray maximum in a transparent SN. 
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Figure 6: Series of the same spectrum as shown in figure 5 (left panel) for different energy bin sizes. From 
top left to bottom right, the bin sizes are 15 keV, 10 keV, 7.5 keV, 5 keV, 2.5 keV, 1 keV. Accounting for the 
high resolution of SPI with its instrumental resolution of «2.3 keV at 850 keV, one single broad Gaussian is 
not representating the observed data any more. Apparently, at least three separate features can be identified 
and followed in time (see below). 

Tracing the 847 keV feature at 10 keV binning in figure 7 (center panel), beginning in epoch (4), 
we identify a broadened line with a velocity spread of (4570 ± 1850) km s” 1 at rest energy, fully 
consistent with the expectations. Towards earlier epochs (3) and (2), we find this line consistently, 
slightly blue-shifted and also broader. Epoch (1), in contrast, shows a red-shifted (f«7000 km s -1 ) 
broad line, together with a narrow blue-shifted feature. Which component corresponds to the bulk 
or peculiar 56 Co emission is arbitrary, and can only be figured out by using smaller time bins, a 
higher spectral resolution and a more sophisticated explosion model: For the 847 keV line it was 
possible to divide the data set to 11 time bins. Using too short time intervals will drown any infor¬ 
mation about the spectral shapes in statistical fluctuations, and makes the detection of the different 
components difficult. The energy band from 780 to 920 keV was then analysed in 2 keV bin¬ 
ning, resolving as narrow as the instrumental line width, for the 11 time bins independently. All 
earlier spectra show more than one spectral feature, and these are all less broad but combine to 
the broad features which are apparent when using coarser energy bins. The last three time bins, 
corresponding to epoch (4), show a consistent picture of a broad 56 Co emission line on top of a 
rather flat spectrum in its vicinity. This reassures the strange behaviour at early times. Before day 
99, «10 keV wide features emerge between 830 and 860 keV, varying in intensity and peak po¬ 
sition with a consistent trend towards lower energies. At «810 keV, another line feature appears 
to rise and fall between days 60 and 90. Another feature at «865 keV which appears to be gone 
in the broad time and energy binning analysis seem to follow a radioactive 56 Co decay. Such de¬ 
tailed discussion may over-interpret our data but we aim to point out a deviation from smoothness. 
We test the hypothesis of the early three epochs being statistical variations of the last, transpar¬ 
ent, epoch by assuming the intrinsic line shape of epoch (4) and fit this template to the earlier 
epochs. The ^ 2 -tests for epoch (2) and (3) reject this hypothesis at the ~ 2 a level. Especially in 
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epoch (1), the discrepancy is very large 
(> 4a), reassuring our finding of asym¬ 
metries in SN2014J. We also have used 
light curves of these three major spectral 
features to test a simple model incorpo¬ 
rating individual 56 Ni contents, kinemat¬ 
ics and occultations. For the components, 
centered at 810, 845 and 865 keV, re¬ 
spectively, we find, consistently, early oc¬ 
culted and later revealed decaying 56 Co 
line emission, and determine optical 
depths at day 1 of 2000, 600, and 100 
for the three components (Maeda, priv. 
comm.). The kinematics derived from 
this three-clump model are a blue-shifted 
«6400 km s -1 (at f«865 keV), a red- 
shifted ss 13000 km s _1 (at «810 keV), 
and a main, spherical symmetric contri¬ 
bution (at ?3 845 keV) to the total 56 Ni 
mass. The relative parts of each feature 
are 15%, 30% and 55%, respectively. We 
do not propose this to be the real mor¬ 
phology of SN2014J, but rather mention 
that this analysis supports our interpreta¬ 
tion of an aspherical 56 Ni distribution due 
to the differences in the respective occul- 
tation of the major features by outer SN 
material. Clumps or co-moving volume 
elements carrying 5b Co may lie along less-occulted lines of sight at specific times. For this reason, 
as the SN expands, different volume elements may thus contribute at different times, as long as 
occulation of the y-rays is significant. 

In a last step we compare the light curve of our brightest 56 Co feature at 847 keV (see figure 8) 
to models, prepared by Lih-Sin The and Adam Burrows [29] specifically for SN2014J, and use 
three different approaches to determine the initially synthesized 56 Ni mass. Details about the pro¬ 
cedure can be found in [38]. The 56 Ni is then determined by a template fit of the models by [29], 
with a global scaling parameter, corresponding to the 56 Ni mass. The most plausible models give 
a 56 Ni mass in the range 0.46-0.59 (±0.06) M 0 , and all imply a substantial fraction of 56 Ni near 
the outskirts of the SN. We note that distinguishing between different models from their fit-quality 
is difficult as the x~ is not varying very much for most of the models. The detailed summary of 
the fitting results can be found in [38]. Taking the average of the three best-fitting models in all 
three approaches yields a 56 Ni mass of (0.49 ±0.09) M 0 . In general, the values range between 
0.4 and 0.8 M 0 . The preference of specific models suggests He on the outside of the SN. These 
(ID) models assume a single broad line to create their time profiles, which of course do not trace 
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Figure 7: SN2014J signal intensity variations for the 
847 keV line (center) and the 1238 keV line (right) as 
seen in the four chosen epochs, and analysed with 10 keV 
energy bins. The 1238 keV fluxes have been scaled by 
1/0.68. Fitted line shape parameters are discussed in the 
text. For the 847 ke V line, in addition, a high-spectral res¬ 
olution analysis is shown at 2 keV energy in width (left). 
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Figure 8: SN2014J signal intensity variations for the 847 keV line in two different time resolutions. The 
4-epoch results are consistent with 11-epoch analysis, both showing an initial rise and late decline of 56 Co 
decay line intensity, with a maximum at 60 to 100 days after the explosion. Shown are also several candidate 
source models from [29], fitted in intensity and thus determining 56 Ni masses for each such model. The 
best-fitting model is shown as a continuous thick line. 


the details in the flux measurements, even for the best-fitted models. If the spectral constraints are 
relaxed, the light curve fits improve, and hence are in good agreement with our proposed irregular 
appearence of the 56 Co lines in SN2014J. We note that in spite of clear detections of the 56 Co y-ray 
lines, discriminating among models is still inadequate, and most models are in general agreement 
with our measurements if 56 Ni is a free parameter. A general model sophistication including dif¬ 
ferent ejecta morphologies, radiation transport and a 3D approach would be in order, to tighten the 
constraints on possible progenitor systems and explosion scenarios. 

5. Conclusions 

We analysed the full set of INTEGRAL/SPI observations available for SN2014J and concen¬ 
trated on the energy ranges around the brightest y-ray line emission expected from the 56 Ni decay 
chain at 158 keV, 812 keV, 847 keV, and 1238 keV. The large instrumental BG of SPI is modelled 
by a detailed spectral analysis which exploits the instrumental spectral response, and accounts for 
continuum and many instrumental y-ray lines arising from cosmic-ray interactions with the in¬ 
strument and spacecraft, in a broader energy range around the lines of interest. In this new BG 
modeling approach, we first identify long-term spectral features, then separately fit their appear¬ 
ances in each individual Ge detector with its response characteristics, and finally determine the 
amplitude variations of the spectral background template on the shorter time scales of our individ¬ 
ual pointings. This method overcomes limitations of statistical precision of BG measurements in 
the fine energy bins where we wish to analyse the source, and uses BG data as measured by the SPI 
detectors themselves, and the nuclear reaction physics extracted from long-term monitoring of BG 
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and detector response (Siegert et al., in prep.). 

Our measurements of 56 Ni decay in SN2014J find the two strongest y-ray lines at 158 and 812 keV 
at day 20 after the explosion. This is surprising not only because of the detection itself but due 
to the kinematics derived. Since most of the 56 Ni y-rays should be absorbed in the SN cloud at 
that time, a major amount of 56 Ni must have been produced outside the SN. The constraints on 
the measured low Doppler-shifts and -velocities indicate an asymmetry in SN2014J. One plausible 
explanation may be a belt of He accreted from a companion star which is exploding, producing 
56 Ni on the surface of the WD, and triggering the SN. 

Moreover, we find emission of decaying 56 Co during the whole observation period of SN2014J at 
847 and 1238 keV, respectively. We analyse separately different epochs of our observations, and 
derive an intensity time profile in the energy regime around the two lines. The brightness evolu¬ 
tion is consistent for both lines, and is overall consistent with expectations from different plausible 
explosion models. We find irregularities in the appearance of the line emission during the rise of 
the measured SN la y-ray light-curve which can not be described as a smooth emergence of the 
Doppler broadened emission lines as they are found in the late, y-ray transparent phase. The de¬ 
rived spectra indicate transient emission features of different widths which is hard to assign clearly 
to 56 Ni or 56 Co. The three-dimensional structure of the SN does not follow the regular shape as 
predicted by one-dimensional models. 

We note that our data cannot discriminate among ID models, although our measurements con¬ 
sistently favour models in which He is on the outskirts. In both measurements we find that the 
explosion SN2014J was not spherically symmetric, dervied from the observed line shapes. These 
results demonstrate the power of y-ray line observations from SNe la, and seed high expectations 
should a closer event occur during the ongoing INTEGRAL mission. 
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